The strength of nuclear shell effects at N = 126 in the r-process region 
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We have investigated nuclear shell effects across the magic number N = f26 in the region of 
the r-process path. Microscopic calculations have been performed using the relativistic Hartree- 
Bogoliubov (RHB) approach within the framework of the Relativistic Mean-Field (RMF) theory for 
isotopic chains of rare-earth nuclei in the r-process region. The Lagrangian model NL-SV1 with the 
inclusion of the vector self-coupling of oj meson has been employed. The RMF results show that 
the shell effects at N = 126 remain strong and exhibit only a slight reduction in the strength in 
going from the r-process path to the neutron drip line. This is in striking contrast to a systematic 
weakening of the shell effects at TV = 82 in the r-process region predicted earlier in the similar 
approach. In comparison the shell effects with microscopic-macroscopic mass formulae show a near 
constancy of shell gaps leading to strong shell effects in the region of r-process path to the drip 
line. A recent analysis of solar-system r-process abundances in a prompt supernova explosion model 
using various mass formulae including the recently introduced mass tables based upon Hartree-Fock- 
Bogoliubov method shows that whilst mass formulae with weak shell effects at TV = 126 give rise to 
a spread and an overproduction of nuclides near the third abundance peak at A ~ 190, mass tables 
with droplet models showing stronger shell effects are able to reproduce the abundance features 
near the third peak appropriately. In comparison, several analyses of the second r-process peak at 
A ~ 130 have required weakened (quenched) shell effects at N — 82. Our predictions in the RMF 
theory with NL-SV1, which exhibit weaker shell effects at N = 82 and correspondingly stronger 
shell effects at N = 126 in the r-process region, support the conjecture that a different nature of the 
shell effects at the magic numbers may be at play in r-process nucleosynthesis of heavy nuclei. 

PACS numbers: 21.10.Dr, 21.30.-x, 21.60.-n, 25.30.+k, 26.50,+x 
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I. INTRODUCTION 

About half of nuclei heavier than Fe are synthe- 
sized in the process of rapid neutron capture (r-process) 
[H S El In environments of high neutron densi- 
ties and high temperatures, extremely neutron-rich nu- 
clei with at least 10-30 neutrons away from the stability 
line are produced. These nuclei arc highly unstable and 
experimentally inaccessible, especially those in the heavy 
mass region. The ensuing nuclei undergo a sequence of 
neutron capture accompanied by a spate of j3~ decays 
thus leading to formation of heavy elements in nature. 
The r-process path passes through the magic numbers 
N = 50, 82 and 126 at different mass values. The syn- 
thesis of nuclei around these magic numbers is reflected 
vividly in the known nuclear abundance peaks around 
A ~ 80, 130 and 190, respectively. 

The shell effects at the magic numbers play a signifi- 
cant role in determining the r-process nuclear abundances 
0. The question whether the shell effects near the r- 
process path are strong or do quench has become crucial 
to understanding the nucleosynthesis of heavy nuclei [(| . 
This question remains open to-date and the existing data 
do not suffice to answer this question. This is due to the 
reason that r-process nuclei are extremely neutron rich 
and are not accessible experimentally. Moreover, on the 
basis of a few nuclei that are known in the extreme re- 
gions, it is not easy to make reliable predictions in the 
farther regions of the period table. Consequently, it is 
proving to be difficult to ascertain the nature of the shell 
effects in the vicinity of the r-process path. A knowledge 



about properties of these nuclei is, therefore, obtained 
from theoretical models. At the same time, new data 
in unknown regions are being obtained experimentally. 
Such data can be of enormous value in defining the na- 
ture of properties of nuclei in the extreme regions. 

In principle, microscopic calculations within a reliable 
model would be attractive for the purpose. The primary 
condition on utility of a microscopic framework should be 
its ability to reproduce features and properties of nuclei 
in the known domains with better accuracy. Calcula- 
tions within a microscopic model for a considerably large 
number of nuclei can be cumbersome. However, with the 
progress in computing speeds, such a task is no more 
beyond one's reach. 

Heretofore, macroscopic-microscopic approaches have 
largely been used to calculate and extrapolate proper- 
ties of nuclei in the inaccessible regions. Most promi- 
nent amongst these is the approach of the Finite-Range 
Droplet Model (FRDM) Q. The mass formula FRDM 
has been obtained on the basis of extensive fits of more 
than a thousand known nuclei across the period table 
including those discovered at the periphery of the pe- 
riodic table in the last decade. R-process calculations 
have been performed using the binding energies (masses) 
and neutron separation energies from the FRDM in con- 
junction with /3-decay properties obtained in the Quasi- 
Particlc Random Phase Approximation (QRPA) Q ■ An- 
other mass formula that has also been employed exten- 
sively is based upon the Extended Thomas-Fermi model 
with Strutinsky Integral (ETF-SI) 9]. Herein, the liquid 
drop (or the smooth) part is provided by the Extended 
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Thomas-Fermi model and the shell corrections are super- 
imposed thereupon using the method of the Strutinsky 
shell correction [lfl| |. 

Employing the results of these two mass formulae ex- 
tensive r-process network chain calculations have been 
undertaken [3- 0- E3| It was concluded that due to strong 
shell effects (gaps) at N = 82 and at N = 126 in the 
region of the r-process path inherent in the mass models 
FRDM and ETF-SI, strong deficiencies (troughs) are ob- 
tained in reproducing solar system r-process abundances 
below the peaks A ~ 130 and A ~ 190 0,0]. A remedial 
measure was suggested by the Hartree-Fock-Bogoliubov 
(HFB) calculations G3 

using the Skyrme force SkP. This 
force is shown to quench (weaken) the strength of the 
shell effects in the r-process region significantly, espe- 
cially at N = 82. This feature of HFB + SkP has proved 
to be useful in filling up the troughs below A ~ 130 and 
A ~ 190 in the r-process abundance curve. In this ap- 
proach, the quenching effect is attributed mainly to a 
large effective mass m* = 1 of the force SkP. Inspired by 
the usefulness of the quenching along the r-process path, 
contrary to the original feature of both the FRDM and 
the ETF-SI, quenching has been introduced in the new 
variant of the mass formula, viz., ETF-SI (Q) [12. The 
mass table ETF-SI (Q) has been shown to be success- 
ful in removing the main deficiencies in the abundance 
curve 13]. These experiments with the mass formulae 
have indicated the need of weaker shell effects along the 
r-process path. A microscopic basis of such a requirement 
at the r-process path, however, needs to be examined. 

In a recent investigation of the r-process nucleosyn- 
thesis of heavy nuclei using mass formulae based upon 
Hartree-Fock-Bogoliubov approach, it has been shown 
[T3 | that weak shell effects in microscopic mass formu- 
lae result in a spread of abundance distribution near 
the A ~ 130 and A ~ 190 peaks. This has the con- 
sequence that large deviations are observed compared 
to the solar-system abundances especially for the peak 
about A ~ 190. In an earlier analysis [l5(, it was also 
shown that in a realistic astrophysical scenario a mass 
model without quenching at N — 126 can fill up defi- 
ciencies (troughs) near A ~ 175 due to freeze-out effects. 
However, this seems to apply to only the third peak in 
the abundance curve. Thus, the role of the shell effects 
at N — 126 in the r-process region is not yet clear. 

In our earlier work 0, 0, 0], we investigated the 
behaviour of the shell effects at the magic number N — 
82 within the Relativistic Hartree-Bogoliubov approach. 
Using the Lagrangian model of the nonlinear quartic cou- 
pling of u> meson in the RMF theory, it was shown that 
microscopic RMF calculations show a weakening of the 
shell effects at N = 82 in going from the stability line 
to the r-process path. It was also shown that in going 
from the r-process path to the neutron drip line, the the 
shell gap diminishes to a vanishingly low value at a given 
isospin, resulting in a complete washout of the shell ef- 
fects In the present work, we have investigated the 
behaviour and evolution of the shell effects at TV = 126 



in the region of the r-process path within the framework 
of the RHB approach. The salient features of the for- 
malism are discussed in Section fllll with emphasis on the 
shell properties of nuclei as inherent to the subject. Re- 
sults of microscopic RMF calculations using two different 
Lagrangian models are presented in Section [3 A com- 
parison of the results is made with the predictions of 
various macroscopic-microscopic mass formulae and in- 
fluence of shell effects on the r-process nucleosynthesis 
is discussed. A discussion is presented of the possible 
consequences of the shell effects on the r-process nucle- 
osynthesis. A summary of the results is presented in the 
last section. 



II. THE SHELL EFFECTS IN NUCLEI 

The shell effects constitute an important feature of nu- 
clei and are known to manifest strongly in terms of the 
magic numbers. This is signified by a conspicuous pres- 
ence of prominent kinks about the major magic numbers 
in two- neutron separation energies (S^n) all along the 
stability line |19| . This is a manifestation of the exis- 
tence of large shell gaps at magic numbers in nuclei. The 
spin-orbit interaction plays a pivotal role in the creation 
of the shell gaps and consequently of the magic numbers. 
Of late, there are indications that with extreme isospin in 
light nuclei, re-adjustment of single-particle levels might 
lead to an emergence of new magic numbers 20] other 
than those which are are hitherto established. 

The spin-orbit interaction and consequently how the 
shell effects behave in the extreme regions would play a 
significant role in carving out shell gaps in nuclei near 
the r-process path. In microscopic approaches such as 
nonrclativistic density-dependent Skyrme theory and the 
relativistic mean-field theory, the spin-orbit interaction is 
determined by data on a few nuclei. Whilst in the former, 
the spin-orbit interaction is added in the theory on an ad 
hoc basis and its strength is adjusted to known spin-orbit 
splittings in a few nuclei, it arises naturally in the RMF 
theory as a consequence of the Dirac-Lorentz structure of 
nucleons. This has shown much usefulness in explaining 
properties that involve shell effects, such as anomalous 
isotope shifts in stable nuclei, especially those associated 
to the Pb chain [2l| . This feature of the shell effects has 
not been possible to attain in the Skyrme theory without 
undertaking significant alterations in the isospin depen- 
dence of the spin-orbit interaction [2^J . On the other 
hand, the intrinsic form of the spin-orbit interaction in 
the RMF theory has been found to be advantageous over 
that in the nonrelativistic approach. The appropriate 
isospin dependence of the spin-orbit interaction |23j has 
been found to be successful in reproducing the anoma- 
lous isotope shifts in Pb nuclei [2l| as well as in Sr and 
Kr isotopes [24| . Consequently, it is expected to have im- 
plications in predicting the shell strength in the extreme 
regions of the r-proce ss p ath. 

The RMF theory [13 [H IH ^ has shown an im- 
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mense potential in being able to describe properties of 
nuclei along the stability line |2^, |3(| and for a large 
number of nuclei beyond the stability line. Most of the 
Lagrangian parameter sets are based upon reproduction 
of binding energies, charge radii and in some cases sur- 
face thickness of a few key nuclei [2f| ■ Various forces are 
obtained in such a way that spin-orbit splitting in some 
key nuclei such as 16 is reproduced reasonably well. It 
should, however, be pointed out that it does not neces- 
sarily ensure that shell gaps or shell effects at the major 
magic numbers are reproduced correctly. 

The shell effects were not considered explicitly in the 
initial developments in the RMF theory of finite nuclei. 
This problem was addressed in Ref. [3JJ, where shell ef- 
fects were investigated in nuclei along the stability line 
with a view to see their influence on nuclei near r-process 
path or on some known "waiting-point" nuclei It 
was shown [3lJ that the otherwise successful RMF forces 
based upon non-linear self-coupling of cr-meson such as 
NL-SH |3(j overestimate the experimental shell gaps in 
nuclei along the stability line. In order to solve this prob- 
lem, the Lagrangian model with the nonlinear scalar cou- 
pling of a meson was extended with the inclusion of the 
nonlinear quartic coupling of the u) meson |3l). Conse- 
quently, shell effects in Ni and Sn isotopes at the stability 
line were reproduced well. 

The shell effects along the stability line (may) have 
repercussions (as it seems to be the case for the RMF 
theory, but not unequivocally for the mass formulae) on 
the shell effects far away from it. The character of the 
shell effects, be it strong or weak vis-a-vis experimen- 
tal data along the line of stability is likely to extrapo- 
late alike (unless major re-adjustments in single-particle 
scheme take place giving rise to unexpected pattern of 
behaviour) in the unknown regions of the periodic ta- 
ble. A test case for this hypothesis was provided by the 
waiting-point nucleus 80 Zn (N = 50) which lies close to 
the r-process path at N = 50. It was shown that 
forces such as NL-SH that overestimate the shell effects 
at the stability line overestimate the shell effects for the 
waiting-point nucleus 80 Zn and also in r-process nuclei 
at N = 82. On the other hand, the force NL-SV1 based 
upon the vector self-coupling of cj-meson, which repro- 
duces the shell effects in nuclei at the stability line, is 
able to reproduce the available data on the waiting-point 
nucleus 80 Zn [l^. A firmer verification of predictability 
of various theories would be provided by future experi- 
mental data in the extreme regions. In the present work, 
we explore how this feature translates for the shell effects 
at N = 126 in the r-process region. 



III. THE RELATIVISTIC MEAN-FIELD 
THEORY 

The RMF approach [23 is based upon the Lagrangian 
density which consists of fields due to various mesons in- 
teracting with nucleons. The mesons include the isoscalar 



scalar cr-meson, the isoscalar vector w-meson and the 
isovector vector p-meson. The details of the formalism 
can be found in Refs. E |H |H H| . 

The RMF Lagrangian that describes the nucleons as 
Dirac spinors moving in the meson fields is given by 

C = 9u<P - g P pf - ^e(l - t 3 )4 - g a a - M N J ip 

+\d^a - U(a) - -fV^ + \mlu^ (1) 
+ ± gi (to^) 2 - \R^ V + - -F pv F^ 

where Mjq is the bare nucleon mass and ip is its Dirac 
spinor. Nucleons interact with a, lj, and p) mesons, with 
the masses being m^, m w and m p and the coupling con- 
stants being g a , g u , g p , respectively. The electromagetic 
interaction is represented by the electromagnetic vector 
field A» . 

The field tensors for the vector mesons are given as 
= d p uj v — dvWp and by similar expressions for the 
p-meson and the photon. For a realistic description of 
nuclear properties a nonlinear self-coupling 

U(a) = \rnla 2 + ±g 2 a 3 + \g^ (2) 

for (7-mesons has been widely used. The non-linear vector 
self-coupling of w-meson as added earlier [2l| m the 
Lagrangian with the non-linear scalar field is represented 
by the coupling constant g4. 

In the lowest order of the quantum field theory, i.e., in 
the mean-field approximation, the nucleons are assumed 
to move independently in the meson fields. The meson 
fields are replaced by their classical expectation values. 
The variational principle leads to the Dirac equation: 

{~iaV + V(r) + /3m*} V,: = e^i (3) 

where V(r) represents the vector potential: 

V(r) = g^ {r) + g p r 3 p (r) + 6(1 ~ Ta) A (r) (4) 

and S(r) is the scalar potential 

S(r) = g a a{v) (5) 

which defines the effective mass as: 

m*(r) =m + S(r) (6) 

The Klein-Gordon equations for the meson fields are 
time-independent inhomogeneous equations with the nu- 
cleon densities as sources. 

{-A + m2}cr( r ) = - 5CT p s (r)- 52 cr 2 (r)- 53( 7 3 (r) 

{-A + m^}w(r) = . 9w p„(r)+.g 4 w 3 (r) 

{-A + m 2 p }p{r) = g p p 3 {r) 
-AA(r) = ep c (r) (7) 
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The stationary state solutions ifa are obtained from 
the coupled system of Dirac and Klein-Gordon equa- 
tions. The ground-state of the nucleus is described by 
a Slater determinant |$ > of single-particle spinors ipi (i 
= 1,2,.. ..A). Solution of the Dirac equation is achieved 
by using the method of oscillator expansion |29j. In the 
RMF approach, the pairing is included within the BCS 
scheme. However, for the case of nuclei in the extreme 
regions of the r-process path and drip lines, the Fermi 
energy is very close to the continuum and many single- 
particle states couple to the continuum. Thus, the BCS 
method of pairing provides a crude approximation of such 
cases. The Relativistic Hartrcc-Bogoliubov (RHB) ap- 
proach based upon quasi-particle scheme provides an ap- 
propriate framework to deal with nuclei of such a nature. 

A. The Relativistic Hartree-Bogoliubov approach 

Nuclei that are known to show strong pairing corre- 
lations are treated appropriately within the framework 
of the RHB approach. Pairing correlations in the neigh- 
bourhood of the Fermi energy in case of nuclei near r- 
process path and drip line become even more important. 
Herein, the RHB approach provides a suitable framework 
to deal with nuclei in the extreme regions. 

It has been shown [34[ that using Green's func- 
tion techniques [j^l a relativistic Hartree-Bogoliubov ap- 
proach can be implemented using the Lagrangian as given 
above. Neglecting retardation effects and the Fock term, 
one obtains relativistic Dirac-Hartree-Bogoliubov (RHB) 
equations 

(-1* -h-)(v) k = E *(v) h > (8) 

where E k are quasiparticle energies and the coefficients 
Uk and Vk are four-dimensional Dirac spinors normalized 
as 

J(U+U k , + V+V k ,)d 3 r = 6 kk ,. (9) 

The average field 

h = a.p + g^uj + 0(M + g a a) - A (10) 

contains the chemical potential A. The meson fields a and 
uj are determined self-consistently from the Klein Gordon 
equations as done in the case of the RMF equations dis- 
cussed above with the scalar density p s = V k V k and 
the baryon density p v — ^2 k V k + V k . The sum on k is 
taken only over the particle states in the no-sea approx- 
imation. The pairing potential A in Eq. JSJ is given by 

A ab = \Y. v «Uc* (n) 

cd 

The RHB equations iJHJ are a set of four coupled integro- 
differential equations for the Dirac spinors U (r) and V(r) 



that are obtained self-consistently. The RHB calculations 
are performed by expanding fcrmionic and bosonic wave- 
functions in 20 oscillator shells. For the pairing channel, 
we use the finite-range Gogny force D1S |3j|. The Gogny 
force is a sum of two Gaussians with finite range. It has 
been shown |3(| that the Gogny force is able to describe 
pairing properties of a large number of finite nuclei in 
the medium and heavy mass regions. Details of the RHB 
theory can be found in Ref. j.3 11 ] 



B. Lagrangian Models 

The Lagrangian model with the nonlinear scalar cou- 
pling of a meson has been the widely used one for finite 
nuclei within the RMF theory. It has been successful in 
reproducing ground-state properties of nuclei at the sta- 
bility line as well as of those far away from it. Here, we 
will consider the successful forces NL-SH [3(j and NL3 
[37j within this Lagrangian model. We will also employ 
the forces NL-SV1 and NL-SV2 [33| (see ref. HH for the 
parameter sets) with the nonlinear vector self-coupling of 
uj meson. As mentioned above, forces NL-SV1 and NL- 
SV2 were constructed with the inclusion of the quartic 
vector coupling of uj meson, in order to solve the prob- 
lem of strong shell effects with Lagr ang ian model with 
nonlinear scalar coupling of a meson . The introduc- 
tion of the non-linear coupling of w-meson also softens 
the equation of state (EOS) of the nuclear matter sig- 
nificantly. This has the consequence that the maximum 
neutron star mass with such an EOS would show a better 
agreement with empirically observed values. A detailed 
discussion of the properties associated with the introduc- 
tion of the nonlinear vector self-coupling of uj meson in 
the RMF theory will be presented elsewhere [3^1 . It will 
be shown [33 that the Lagrangian parameter set NL-SV1 
is also able to improve upon the ground-state properties 
such as binding energies, charge radii and isotopes shifts 
of nuclei along the stability line and far away from it as 
compared to those with NL-SH and NL3. 



IV. DETAILS OF THE CALCULATIONS 

In RHB calculations, wavefunctions are expanded into 
an harmonic-oscillator basis to solve the Dirac and the 
Klein-Gordan equations [2^| . For both the fermionic and 
bosonic fields a basis of 20 oscillator shells has been used. 
The pairing has been taken in the Bogoliubov approach. 
For, nuclei a few neutrons below and a few neutrons 
above a magic number are usually spherical, RHB calcu- 
lations have been performed for a spherical configuration. 
For a comparative study of the shell effects, we have per- 
formed RHB calculations with the two Lagrangian mod- 
els as discussed above. However, our focus is on investi- 
gation of the potential and predictive power of the new 
Lagrangian model with the vector self-coupling of uj me- 
son vis-a-vis the scalar self-coupling of a meson. 
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With a view to investigate as to how the shell effects 
evolve in going from the region of the stability line to- 
wards the r-process path and ultimately to the neutron 
drip line, we have selected even-even nuclei from the iso- 
topic chains of Hf [Z = 72) down to Xe (Z = 54) across 
the neutron magic number TV = 126. For our focus is on 
the behaviour of the shell gap, nuclei relevant to the dis- 
cussion are those with N = 124, N = 126 and N = 128, 
so as to be able to calculate Sz n values across the magic 
number N = 126. 



V. RESULTS AND DISCUSSION 
A. Shell effects along the stability line 

First, we examine the known shell effects at N = 126 at 
the stability line. The experimental shell gap at N = 126 
is known only in a very few isotopic chains in this region. 
Here, the doubly closed nucleus 208 Pb in the Pb chain 
provides a cardinal point in the study of the shell effects 
at N = 126. Interestingly, though the nuclear landscape 
of the periodic table has been extended significantly due 
to sustained experimental efforts in the last decade to- 
wards synthesizing nuclei far beyond the stability line 
in the laboratory, the heaviest known Pb nucleus has 
reached only to 214 Pb. 

We begin with the premise that the two-neutron sepa- 
ration energy at the magic number provides a reasonably 
good indicator of the shell gap. Therefore, we calculate 
the shell gap at the magic number as defined conveniently 
by 

A s = S 2n (Z,No)-S2n(Z,N + 2), (12) 

where S2u(Z,Nq) denotes the 2-neutron separation en- 
ergy of the nucleus (Z, Nq) with a magic neutron number 
N . 

The shell gap in Pb nuclei obtained from the RMF 
approach using various Lagrangian sets is shown in 
Fig. 1(a) and (b). The corresponding gap from vari- 
ous macroscopic-microscopic mass formulae is shown in 
Fig. 1(c). The results are compared with the experimen- 
tally known data. The S2n values calculated in the RMF 
theory with the forces NL-SH and NL3 within Lagrangian 
model of the nonlinear scalar coupling of a meson are 
shown in Fig. 1(a). The difference between the data 
points at N = 126 and N = 128, as indicated by Eg. ITU 
manifests the shell gap at the magic number N = 126. 
A comparison with the experimental data points shows 
that the shell gap from the force NL-SH underestimates 
the experimental gap by ~ 0.5 MeV. On the other hand, 
the recent force NL3 shows a good agreement with the 
data. Here, we do not show the binding energy of 208 Pb 
itself that is overestimated both by NL-SH and NL3 by 
about 2-3 MeV. However, the difference in the S^n values 
of the neighboring nuclei turns out to be satisfactory in 
both the cases. It is interesting to note that although 
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FIG. 1: The two-neutron separation energies Sin for Pb iso- 
topes as obtained from (a) the forces with the Lagrangian 
model of the nonlinear scalar self-coupling of a meson, (b) 
from the forces which include the nonlinear quartic coupling of 
ui meson, and (c) from the mass formulae FRDM 0, ETF-SI 
[9|l and the recently obtained results from mass tables HFB-2 
|39|. The experimental data is shown by the solid circles. 



NL-SH underestimates the shell gap at N — 126 slightly, 
it is known to exhibit generally stronger shell effects for 
nuclei at and beyond the stability line than those from 
NL3. This implies that a focus on a single or a few data 
points may not be deterministic as far as extrapolations 
(predictions) for nuclei beyond the stability line are con- 
cerned. Consequently, the shell gap at N = 126 in the 
Pb nucleus may not provide a successful conjecture as to 
whether the shell effects with either of the forces shall 
remain strong or weak in the domain that is far beyond 
the stability line. 

The Sm values obtained from the forces NL-SV1 and 
NL-SV2 with the Lagrangian model with the nonlinear 
quartic coupling of to meson are shown in Fig. 1(b). A 
comparison with the experimental data shows that NL- 
SV1 reproduces the shell gap at N = 126 well, whereas 
it is underestimated by ~ 0.5 MeV by the force NL-SV2. 
Again, a seemingly paradoxical situation arises here. As 
in the case of NL-SH, exhibiting generally stronger shell 
effects and yet underpredicting the shell gap at N = 126 
as shown in Fig. 1(a), the force NL-SV2 has also been 
shown to exhibit shell effects slightly stronger than those 
with NL-SV1 in the r-process region at N — 82 [l6| . 
This shows that reliance on a few data points may not 
be useful in predicting behaviour in the extreme regions. 
This will be shown in the latter parts of this paper, where 
we will discuss the shell effects in the region of the r- 
process path and the neutron drip line. 

It is equally interesting to see as to how various mass 
formulae predict the shell gap at N = 126 in Pb nuclei 
at the stability line. In Fig. 1(c), the data points from 
the mass formula FRDM and ETF-SI are shown. It may 
be remarked that various mass formulae including FRDM 
and ETF-SI have been obtained with a view to reproduc- 
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FIG. 2: The two-neutron separation energies S2n for the isotopic chains from Hf [Z = 72) to Xe (Z = 54) in going towards the 
r-process nuclei and in approaching the neutron drip line, obtained from the forces with the Lagrangian model of the nonlinear 
scalar self-coupling of a meson (a) NL-SH [30| and (b) NL3 |37jl . The 5*2n values from the forces which include the nonlinear 
quartic coupling of co meson (3lll33|| are also shown in (c) with NL-SV2 and in (d) with NL-SV1. 



ing experimental data on more than a thousand nuclei. 
Generally, these mass tables have achieved a great suc- 
cess in reproducing a large set of experimental database 
through exhaustive fits over the periodic table than pos- 
sibly a microscopic theory could ever do. However, as 
pointed out in the literature 0, l39| , discrepancies at the 
magic numbers do remain a significant drawback. 



The FRDM indicates a shell gap that is ~ 1.5 MeV 
smaller than the experimental one. On the other hand, 
ETF-SI underestimates the shell gap only by about 0.5 
MeV. The undervaluation of the shell gap with FRDM 
and ETF-SI at N = 126 along the stability line is to 
be contrasted to the stronger shell effects due to these 
mass formulae when extrapolated in the extreme regions 
of the r-process path both at N — 82 and N — 126, than 
are suggested for a successful reproduction of r-process 
abundances. Notwithstanding the need of weaker shell 
effects along the r-process path, a new mass table based 
upon the Skyrme Hartree-Fock-Bogoliubov (HFB) ap- 
proach has recently been produced [3!j . The data points 
from the mass table HFB-2 as shown in Fig. 1(c) underes- 
timate the shell gap at N = 126 significantly. It is worth 
pointing out that the mass formula HFB-2 [3!j seems to 
have achieved a similar quality of fit across much of the 
periodic table. It is comparable to ETF-SI and FRDM, 
albeit with shell gaps that are predicted to be small in 
the r-process region, as we will see in Section IV CI 



B. Shell effects near the r-process path - RMF 
theory 

The shell effects at N = 126 for r-process nuclei play a 
crucial role in determining r-process abundances around 
the peak at A ~ 190 For practical purposes r-process 
nuclei are defined to be those with S n ~ 2 — 4 MeV. 
Thus, the r-process path is not strictly well defined and it 
does vary from model to model. However, it is generally 
accepted that nuclei with Z ~ 64 — 69 near N — 126 fall 
along the r-process path. 

The results of RHB calculations for the forces NL-SH 
and NL3 with the Lagrangian model of nonlinear a cou- 
pling are shown in the upper panels (a) and (b) of Fig. 2, 
respectively. It is seen clearly that the shell gap that 
is represented by the difference between the Sm values 
at N = 126 and at N — 128 shows a gradual decrease 
with an increase in isospin i.e., in going from the element 
Hf (Z = 72) that is slightly above the r-process path to 
Xe (Z = 54) that is near the drip line. This behaviour 
is similar for both NL-SH and NL3. In absolute terms 
the shell gaps for NL-SH are slightly larger than those 
of NL3. However, in both the cases, the shell gap does 
not show a significant decrease in approaching the drip 
line, as is probably expected from a comparison with the 
correspondin g b ehaviour of the shell gap at N — 82 at 
the drip line |lq . 

The S2n values obtained from the forces NL-SV1 and 
NL-SV2 with the inclusion of quartic coupling of uj meson 
are shown in the lower panels (c) and (d) of Fig. 2, respec- 
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FIG. 3: The two-neutron separation energies 5*2n for the isotopic chains from Hf (Z = 72) to Xe {Z = 54) in going towards 
the r-process nuclei and in approaching the drip line, as predicted by the mass models (a) FRDM Q and (b) ETF-SI @. The 
S2n values from the mass formula ETF-SI (Q) O that includes a quenching of the shell effects superimposed on ETF-SI arc 
shown in (c). The results from the mass formula HFB-2 [3Sj are shown in (d). 



tively. For both the forces NL-SV1 and NL-SV2, the shell 
gaps show a similar gradual decrease in going towards the 
r-process path and the drip line. Qualitatively, the be- 
haviour of both the Lagrangian models as portrayed in 
the upper and the lower panels, respectively, is very simi- 
lar. However, quantitatively, the shell gaps with NL-SV1 
are smaller than both with NL-SH and NL3. The force 
NL-SV2 provides slightly larger shell gaps than those NL- 
SV1. 

All the parameters sets of both the Lagrangian mod- 
els exhibit a slight reduction in the shell strength in the 
region of the r-process path. The difference lies only in 
the degree by which the shell effects are reduced in going 
from the r-process path to the drip line. We will dis- 
cuss the comparative behaviour of the shell effects with 
various RMF models in Section IV Dl 



Shell effects near the r-process path 
models 



In the absence of and in essence rather infeasibility 
at present of constructing a mass table based purely on 
microscopic calculations, masses from various mass tables 
based upon macroscopic-microscopic approach are used 
in r-process calculations. Here, we present the results of 
the two most elaborate mass formulae, the FRDM and 
the ETF-SI. 

We show the S2n values obtained from the FRDM and 
ETF-SI in the upper panels (a) and (b) of Fig. 3, re- 



spectively. In contrast to the microscopic calculations 
of Fig. 2, the shell gaps with FRDM and ETF-SI show 
nearly constant values in going from Z = 72 to Z = 54. 
On the other hand, a slight increase in the shell gap is 
visible with FRDM in going towards Xe. However, as the 
S 2n values from FRDM for the nuclides with N = 128 
become negative below Nd (Z = 60), the apparent in- 
crease in the value of the shell gap for nuclei below Nd 
can be therefore be discounted. 

A constant shell gap is also displayed by ETF-SI as 
shown in Fig. 3(b). The magnitude of the shell gap with 
ETF-SI, is, however, larger than that with FRDM by ~ 
0.8 MeV. On the lower end, i.e, for nuclides with N = 
128, the ETF-SI also shows negative S^n values below 
Nd (Z = 60). Thus, the drip line is reached at N = 128 
for the elements below Nd, the behaviour very similar to 
that of FRDM. Moreover, with ETF-SI the S 2n values 
for nuclides with the magic number (N = 126) are ~ 1 
MeV larger than those of FRDM. The similar behaviour 
of the shell effects with FRDM and ETF-SI and arrival 
of the drip line at a similar location is not surprising, for 
the shell corrections superimposed on the smooth part in 
the two models are based upon the same prescription of 
the Strutinsky shell correction 

Predictions from the mass formulae FRDM and ETF- 
SI in the extreme regions of the r-process path have 
been used extensively for network chain calculations of r- 
process nuclear abundances. Results of calculations have 
shown that due to strong shell effects that are preva- 
lent with FRDM and ETF-SI along the r-process path, 
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FIG. 4: The S^n values from NL-SV1 are compared to those from (a) NL-SH and (b) NL3. The shell gap region is bounded 
by the two vertical lines at N — 126 and N = 128. 



there is a significant deficiency (troughs) in the r-process 
abundances below the A ~ 130 and A ~ 190 peaks Q. 
Inspired by the Hartree-Fock-Bogoliubov (HFB) calcu- 
lations with the Skyrme force SkP |ll|. quenching was 
introduced in the ETF-SI mass formula, thus creating a 
new mass table ETF-SI (Q) [12|. The resulting r-process 
calculations p| with ETF-SI (Q) have been able to fill 
up the deficiencies below the two peaks and results seem 
to be promising. This has put up a requirement (indi- 
rect) of a weakenin g of the shell structure for nuclei near 
the r-process path | l3t |40| . This feature that is clearly 
absent in the original mass formulae FRDM and ETF-SI 
has been introduced in the ETF-SI (Q) rather artificially. 
The ETF-SI (Q) results in Fig. 3(c) show that the shell 
gaps remain the same as that with ETF-SI in going from 
Hf [Z = 72) to Er [Z = 68). However, the shell gap 
starts decreasing rapidly below Dy (Z — 66) in going 
towards Xe due to the quenching introduced therein. A 
microscopic basis of the aforesaid quenching introduced 
in the mass formula ETF-SI (Q) is yet to be established. 

Motivated by the quenching present in the HFB calcu- 
lations with SkP, attempts have been made to introduce 
the HFB approach in some mass tables. We show in 
Fig. 3 (d) the results taken from the recently developed 
mass table HFB-2 39] within the Skyrme Ansatz. The 
shell gap with HFB-2 is reduced for all the isotopic chains 
as compared to FRDM and ETF-SI. A reduction in the 
shell gap with HFB-2 was also seen for Pb isotopes in 
Fig. 1(c). Though the shell gap is reduced vis-a-vis other 
mass formulae, it, however, remains constant in going 
from Hf [Z = 72) to Xe (Z = 54). There are no indica- 
tions of an additional reduction in the shell strength near 



the drip line as compared to the region of the r-process 
path. Thus, the constancy of the shell gap in going from 
the region of the r-process path to the drip line seems 
to be the salient feature of the macroscopic-microscopic 
mass tables presented in Fig. 3. Here, the only exception 
is ETF-SI (Q), wherein the quenching was introduced by 
force. In comparison, the RMF results show a slight re- 
duction in the shell strength in going from the r-process 
path to the drip line. 

An earlier version of the HFB mass table, i.e., HFB-1 
was constructed by replacing the BCS pairing by the Bo- 
goliubov pairing scheme in Ref. |4l) . In this work, it was 
shown that the behaviour of shell gaps far away from 
the stability line does not depend much upon whether 
the BCS pairing or the Bogoliubov pairing is used. Ac- 
cordingly, the shell gaps at N = 126 with HFB-1 were 
found to be comparable to those with FRDM and that 
an introduction of Bogoliubov pairing did not result in a 
quenching of the shell effects [4l| . 



D. A comparative analysis of the shell effects 

In order to understand the comparative behaviour of 
the shell effects first in the RMF theory, we show in Fig. 4 
the S^n values across the magic number N — 126 from 
the force NL-SV1 and compare these with the other forces 
such as NL-SH and NL3. It may be reminded that al- 
though the force NL-SH has been found to be successful 
in reproducing binding energies, charge radii and defor- 
mation properties of a large number of nuclei far away 
from the stability line, including the anomalous isotope 
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shifts in Pb nuclei, the shell effects with NL-SH were 
found to be stronger as compared to the experimental 
data [Ul • I n comparison, the force NL-SV1, having been 
able to describe the shell effects in Ni and Sn isotopes 
at the stability line, was also shown to be successful in 
reproducing the available data on the shell effects at the 
waiting-point nucleus 80 Zn |l6j |. In view of this, we treat 
the force NL-SV1 as our benchmark. 

In Fig. 4(a), we compare the results of NL-SV1 to those 
from NL-SH. The variation in the shell gap is illustrated 
by the changing slopes of the curves between the vertical 
bars at N = 126 and N = 128. A look at the differ- 
ence in the S-zn values and at the corresponding slope 
of the curve for Hf (Z = 72) shows that the shell gap 
with NL-SH is bigger than with NL-SV1. As one pro- 
gresses towards the r-process nuclei such as Er, Dy, Gd 
and Sm, this difference in the shell gap between NL-SH 
and NL-SV1 increases with an increase in the isospin. 
One sees that the shell gaps with NL-SH remain stronger 
even in nuclei near the drip line such as Xe. In compar- 
ison, NL-SV1 shows a faster decrease in the shell gap in 
going from r-process nuclei to the drip line, a feature that 
has been called for for reproduction of r-process nuclear 
abundances. 

The results from NL-SV1 are compared with those 
from the force NL3 in Fig. 4(b). A comparison between 
the shell gaps from the two forces shows that beginning 
with Hf, the shell gap with NL3 is slightly larger than 
that with NL-SV1. This difference, however, increases 
slowly when one moves from Hf to Xe. Thus, the shell 
effects with NL3 are slightly stronger than those with 
NL-SV1. Additionally, the S-zn values for nuclides with 
N = 124 are - 0.5 MeV higher with NL3 than NL-SV1, 
especially in the region of r-process nuclei. This feature 
is similar to that of NL-SH vis-a-vis NL-SV1 as shown 
in Fig. 4(a). We compare in Fig. 5 the shell gaps from 
NL-SV1 to those from the mass models (a) ETF-SI and 
from its variant (b) ETF-SI (Q) that embeds a quench- 
ing near the r-process region. For Hf [Z = 72), the shell 
gap with NL-SV1 is ~ 1 MeV larger than that of ETF- 
SI. However, as one proceeds to the r-process nuclei and 
towards the drip line, this difference between NL-SV1 
and ETF-SI decreases and then it reverses. For drip line 
nuclei near Xe (Z — 54), the NL-SV1 shell gap is then 
smaller than that with ETF-SI. Whereas the shell gap 
with ETF-SI hardly shows any change in going from Hf 
to Xe, the shell gap with NL-SV1 does show a consistent 
decrease in going from the r-process to the drip line. 

We also compare the NL-SV1 predictions to those of 
ETF-SI (Q) in Fig. 5(b). From the nuclei of Hf [Z = 72) 
to about Sm [Z = 62), there is not much difference be- 
tween the results of ETF-SI and ETF-SI (Q). Therefore, 
for these nuclei, a comparison of NL-SV1 shell gaps with 
ETF-SI (Q) ones is similar to that with ETF-SI as shown 
in Fig. 5(b). In the region of r-process nuclei Z = 64 — 68, 
the shell gaps between NL-SV1 and ETF-SI (Q) are sim- 
ilar. However, due to an extra quenching added in ETF- 
SI (Q), differences in the shell gaps of the two approaches 



begin appearing below Z=62. The shell gap with ETF- 
SI (Q) becomes especially small for nuclei near the drip 
line. The impact of the variations in the shell effects 
in nuclei can not be visualized without comprehensive 
r-process calculations. In the case of ETF-SI (Q), such 
calculations do exist and have shown promising results. 
On the other hand, r-process calculations using the re- 
sults from NL-SV1 are being planned currently. 



E. Models with strong shell effects 

The strength of the shell effects at the major magic 
numbers has been a point of numerous discussions in 
respect of nuclear abundances 0, IE El • The focus has 
mostly been on the possibilities and potential capabili- 
ties of various mass formulae of macroscopic-microscopic 
origin. This is evidently due to the fact that mass formu- 
lae have been able to produce large-scale nuclear bind- 
ing energies and other data relevant for use in network 
chain calculations. The microscopic theories have not 
yet enjoyed this privilege with the exception of isolated 
case(s) where nuclear masses for the purpose have been 
calculated meaningfully. Here, we wish to take stock of 
the situation on the shell effects of nuclei within various 
frameworks. 

We compare results of the mass formulae FRDM and 
ETF-SI, both of which show strong shell effects along 
with those from the RMF force NL-SH also exhibits 
strong shell effects. The S2n values from FRDM and 
NL-SH are shown across the shell closure N — 126 in 
Fig. 6(a). The parallel lines for S 2 n with FRDM from 
Hf to Xe in the the region of the shell gap bounded by 
the vertical lines in Fig. 6(a) indicate the constancy of 
the shell gap with FRDM as discussed above. The shell 
gap as indicated by the £2™ values from NL-SH and the 
ensuing large slope of the curve between N — 126 and 
N = 128 demonstrates overly strong shell effects with the 
force NL-SH. In comparison, the shell gap with FRDM 
is smaller than that with NL-SH by about 2 MeV for Hf 
nucleus. This difference between NL-SH and FRDM per- 
sists even for r-process nuclei in the region of Z — 62 — 66. 
Thus, the shell effects at N = 126 with NL-SH are signif- 
icantly stronger than those with FRDM. However, due 
to natural decrease in the shell gap with an increase in 
isospin in the RMF theory, the shell gap with NL-SH 
does become comparable to FRDM for nuclei near the 
drip line. 

The indication that the shell strength with NL-SH is 
strong appeared in Ref . [42( , where the shell effects with 
NL-SH were studied across N — 82. In that work, it 
was shown that the shell effects at TV — 82 with NL-SH 
were as strong as those with FRDM, but not stronger 
than FRDM. It was also shown that there was a remark- 
able agreement between the ground-state properties of 
the isotopic chains of Zr from 112 Zr through to 130 Zr in 
the two approaches. Nuclei in this region exhibited not 
only similar values of quadrupole deformation, but also 
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FIG. 5: The S2n values from NL-SV1 are compared to those from the mass formula (a) ETF-SI and with the quenched mass 
formula (b) ETF-SI (Q). The shell gap region is bounded by the two vertical lines at N = 126 and N = 128. 




FIG. 6: The S2n values from the FRDM are compared to those (a) from the force with strong shell effects NL-SH and (b) the 
mass formula with strong shell effects ETF-SI. 



similar shape transitions across the region with FRDM 
and NL-SH 0. Tnis to a surmise that the shell ef- 
fects and shell structure with NL-SH are very similar to 
those of FRDM. However, this does not seem to be the 
case in the region of N = 126. 

Taking into consideration the results of network chain 
calculations performed thus far, we believe that r-process 
network chain calculations with the force NL-SH or any 



similar microscopic force exhibiting strong shell effects 
are not expected to be successful, at least for the third 
peak at A ~ 190. 

We compare the shell effects with FRDM to those with 
ETF-SI in Fig. 6(b). The shell gaps with ETF-SI for 
nuclei in the vicinity of Hf are more than 1 MeV larger 
than those with FRDM. The strong shell gaps with ETF- 
SI are maintained across the r-process region and a near 
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FIG. 7: The shell gaps A s (Eq. at N = 126 as obtained 
from the RMF forces (a) NL-SH, NL-SV1 and NL-SV2 and 
compared with those from ETF-SI (Q). (b) The shell gaps 
from mass formulae ETF-SI, FRDM and HFB-2 are shown. 
A comparison is also made with the ETF-SI (Q). 



constancy of the shell gap at N = 126 both with FRDM 
and ETF-SI is seen clearly. The difference in the shell 
gap between ETF-SI and FRDM, however, decreases in 
going from the r-process path nuclei towards the drip 
line nuclei. The fact remains that with ETF-SI the shell 
effects are stronger than those with FRDM. 

Using macroscopic-microscopic mass formulae FRDM 
and ETF-SI, it was shown |1 H El that due to the 
stronger shell gaps in ETF-SI, r-process network chain 
calculations lead to a greater deficiency (trough) in the 
r-process nuclear abundances below A ~ 190 peak than 
with FRDM. Thus, the earlier conclusion of much of these 
analyses has been that stronger shell effects at N = 126 
do not seem to be conducive to reproducing nuclear abun- 
dances. However, it was shown in a later analysis ^3 
that in a "realistic" astrophysical scenario, there is no 
stringent need for a quenching of the N — 126 shell ef- 
fects. The trough that appeared near A ~ 175 in ear- 
lier analyses could be filled due to freeze-out effects even 
by using a mass model without quenching. This raises 
the possibility that under appropriate conditions, mass 
formulae without a quenching of the shell strength at 
N = 126 can be used successfully. 



The N = 126 shell gap in nuclei: RMF versus 
mass formulae 



The status of JV = 126 shell gaps As as defined in 
Eq-Elis summarized in Fig. 7 for the RMF theory and 
for various mass formulae. The shell gap at N — 126 as 
obtained from the RMF calculations with NL-SH, NL- 
SV1 and NL-SV2 are shown in Fig. 7(a). All the micro- 
scopic calculations show a decreasing shell gap in going 



through the r-process region. This reduction is, however, 
a little stronger as one moves towards the drip line. 

A comparison shows that NL-SH exhibits shell gaps 
which are consistently larger than those with NL-SV1 
and NL-SV2. This characterizes the stronger shell effects 
of NL-SH vis-a-vis other forces. The shell gaps with NL- 
SV2 are slightly larger than those with NL-SV1. This 
was also shown to be the case for shell gaps at N = 82 
with NL-SV2 Qjj. 

We compare the shell gaps with the RMF forces to 
those from the quenched mass formula ETF-SI (Q) in 
Fig. 7(a). Evidently, the shell gaps with ETF-SI (Q) are 
much smaller than those with NL-SH. However, these are 
also smaller than those of the NL-SV1 and NL-SV2. The 
effect of the added quenching in ETF-SI (Q) is apparent 
for nuclei below Z = 64, where the shell gap is reduced 
significantly as compared to the nearly constant values 
maintained in ETF-SI (see Fig. 7(b)). For nuclei below 
Nd (Z = 60), however, the quenching in the ETF-SI (Q) 
is much stronger than the weakening of the shell effects 
predicted by NL-SV1. 

The shell gaps from various mass formulae are com- 
pared in Fig. 7(b). Both ETF-SI and FRDM exhibit shell 
gaps which remain nearly constant in going from the re- 
gion of the r-process path towards the drip line. Compar- 
atively, ETF-SI shows shell effects that arc stronger than 
those with FRDM in much of the region shown in the fig- 
ure. As shown in earlier calculations ^|j the strength of 
the shell effects with ETF-SI was found to produce much 
larger deficiency below A ~ 190 peak in the r-process 
abundances. 

The shell gaps with the new mass formula HFB-2 
are compared with those from ETF-SI and FRDM in 
Fig. 7(b). The shell gaps with HFB-2 are nearly con- 
stant as shown also by the other mass formulae. These 
are, however, systematically smaller than those of ETF- 
SI and FRDM. There is a tendency of only a slight de- 
crease in the shell gap in going towards the drip line 
nuclei. It is interesting to note that the addition of a 
Bogoliubov based pairing in the mass formula has not 
been found to be sufficient to suppress further the shell 
gaps below the r-process region. A comparison of the 
shell gaps from ETF-SI, FRDM and HFB-2 with those 
from ETF-SI (Q) shows that the shell gaps from ETF- 
SI (Q) are in striking contrast to all the other mass tables. 
However, as mentioned earlier, this striking difference has 
been brought about by the introduction of a quenching 
based upon results of HFB+SkP calculations. This ex- 
tra weakening of the shell gaps in the drip line region 
is not shown by any of the widely used mass formulae. 
It has been reported in some calculations p"2l E| that 
a weakening of the shell effects in the r-process region 
to the drip line is required for reproducing the r-process 
abundances around the peaks at A ~ 130 and A ~ 190. 
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FIG. 8: The S^n values from NL-SV1 are shown for nuclei 
across the shell closure (a) N — 82 and (b) N = 126 in the 
region of the r-process path. The evolution of the shell gaps 
at N = 126 is compared to that at N = 82. 



G. The shell N = 126 versus N = 82 

We discussed the evolution of the shell gap at TV = 82 
near the r-process path in the RMF theory in Ref. [ly in 
detail. The shell gaps at N — 82 were found to weaken 
in the region of the r-process path. In going to the drip 
line, the shell effects at TV = 82 showed a substantial 
reduction in the shell strength, eventually leading to a 
complete disappearance of the shell gap for the drip line 
nuclei. 

We show in Fig. 8 a comparison of the shell effects at 
N = 126 as obtained in the present work with those at 
N = 82 ^(| in the region of the r-process path to the 
drip line. In order to visualize the shell gap at N = 82, 
we show the S^n values obtained with NL-SV1 across 
the magic number for nuclides of Sn (Z = 50) down to 
the drip line nuclides of Kr (Z = 36). The shell gap 
at N = 82 shows a strong decrease in going from Cd 
(Z = 48) to Zr (Z = 40) in the region of the r-process 
path as mentioned above. This reduction in the shell gap 
at N — 82 is significantly faster than the corresponding 
shell gap at N = 126 (Fig. 8(b)). In the latter case, the 
shell gap shows a much slower decrease in the region of 
the r-process nuclei from Er (Z = 68) to Nd (Z — 60). 
The N — 12& shell strength exhibits a resilience to the 
change in isospin. 



H. The shell effects and r-process nucleosynthesis 

In the broader context of the shell effects and their 
implications on the r-process nucleosynthesis, it is per- 
tinent to discuss recent r-process calculations which em- 
ploy newly developed mass formulae HFB-2 and HFB-7 
|l4j . In this work, effect of various mass formulae on r- 
process nucleosynthesis has been studied using astrophys- 



ical model of prompt supernova explosion from a collaps- 
ing O-Ne-Mg core [44| . In the mass formulae HFB-2 and 
HFB-7 which have been used, pairing has been taken into 
account by the Bogoliubov method in a Skyrme density- 
functional approach. R-process calculations in this study 
have also been carried out using the nuclear masses from 
FRDM 7] and from the old droplet mass (DM) model 

It has been shown 01 that due to the shell effects 
which are significantly reduced (weak) with HFB-2 and 
HFB-7, the abundance curve in the third peak at A ~ 190 
is spread leading to a shift of the peak and consequently 
the valley at A ~ 183 is shifted significantly to lower 
masses. The results show that by the use of HFB-2 and 
HFB-7 masses, there are large deviations in the third 
peak as compared to the solar-system r-process abun- 
dances. Thus, due to the weakness of the shell effects at 
N = 126 near the r-process in the mass models HFB-2 
and HFB-7, there is a significant overproduction of nuclei 
to the left of the third peak in the solar pattern. It is 
interesting to note that on using masses due to FRDM 
and DM, which are characterized by stronger shell effects 
at N — 126, it has been shown [l4| that the abundance 
curve gives rise to a sharp r-process peak at A ~ 190, 
in better agreement with the solar pattern in the re- 
gion. These results suggest that in the prompt supernova 
model considered in Ref. Q and conditions applicable 
therein, mass models/microscopic theories with stronger 
shell effects at N — 126 would reproduce the main fea- 
tures of the solar r-process abundance curve around the 
third peak reasonably well. It may be recalled that rather 
similar conclusion was reached in Ref. |l5j . wherein 
it was shown that a mass model without quenching at 
N = 126 can fill up the trough at A ~ 175 and repro- 
duce the abundance curve near the third peak due to 
freeze-out effects. 

In view of the different behaviour of the shell effects 
at N = 82 and N = 126 as predicted by NL-SV1 in 
the RMF theory as shown in Fig. 8, it is pertinent to 
discuss the results of the r-process calculations of Ref. 
[lH around the second peak at A ~ 130. The re- 
sults have shown that droplet mass models FDRM and 
DM with strong shell effects also at N = 82 produce 
troughs (underproduction) at A ~ 115 and A ~ 140 
in the abundance curve. This shows that in contrast 
with N = 126, stronger shell effects at N = 82 are not 
desirable for reproducing the solar abundance pattern. 
On the other hand, with HFB-2 and HFB-7 deficiencies 
in the r-process production below and above the second 
peak at A ~ 130 are significantly remedied especially 
with HFB-2. However, the weak shell effects inherent 
in HFB-2 and HFB-7 have the consequence that abun- 
dances around A ~ 130 are also spread out as opposed 
to the solar system r-process abundances. 

The results of Ref. ^4| indicate that within the model 
employed in this work, the two peaks at A ~ 130 and 
A ~ 190 require a different nature of the shell effects. For 
a successful reproduction of abundances near A ~ 130, 
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shell effects at N = 82 in the r-process region should 
not be strong, whereas the third peak seems to require 
moderate to strong shell effects at N = 126 in contrast 
to N = 82. Thus, in the model considered in Ref. [T^| . 
the analysis of the second and the third peaks in the 
solar-system r-process abundance curve suggests a differ- 
ent nature of the shell effects at N — 82 and at N = 126. 

Microscopic RHB calculations with the Lagrangian 
model NL-SV1 show two very different features for the 
shell effects at N = 82 and N = 126, as depicted in 
Fig. 8. These features seem to be consistent with the 
picture that emerges from the results of Ref. ^4|. It 
remains to be seen in future network chain calculations 
with masses obtained from microscopic calculations with 
NL-SV1 whether the two different features of the shell 
effects exhibited by NL-SV1 would suffice towards repro- 
ducing the r-process abundances near the second and the 
third peak. 

VI. CONCLUSIONS 

We have investigated the shell effects in nuclei at the 
magic neutron number N — 126 in the region of the 
r-process path using the relativistic Hartree-Bogoliubov 
approach within the relativistic mean-field theory. Two 
Lagrangian models, one with the nonlinear scalar cou- 
pling of <7 meson and another one that includes a non- 
linear vector coupling of u) meson have been considered. 
A comparison of the shell effects in the RMF theory has 
been made with the predictions of various mass formu- 
lae. It is shown that the predictions of RMF calculations 
exhibit a slight reduction of the shell gap in going from 
the r-process path to the neutron drip line irrespective 
of the shell strength exhibited by a Lagrangian parame- 
ter set along the stability line. This is slightly different 
from a near constancy of the shell gaps demonstrated by 
major mass formulae in the region of the r-process path 
and the drip line. Consequences of the shell strength on 



the r-process nucleosynthesis have been discussed. 

It is shown that the Lagrangian force NL-SV1 with the 
vector self-coupling of co meson, which reproduces the 
shell gaps along the stability line, shows that the shell 
effects at N = 126 exhibit only a marginal reduction in 
the shell strength in going from the r-process path to the 
drip line. Consequenty, the shell effects retain a strong 
character along the r-process path at N = 126. This is in 
striking contrast to the earlier results ^(| with NL-SV1 
that in the RMF theory shell effects at N = 82 exhibit 
a significant weakening of the strength in going from the 
r-process path to the neutron drip line. This shows that 
different magic numbers may exhibit a different nature 
of the shell effects in the extreme regions of the periodic 
table. 

Analysis of the results of a recent r-process calculations 
[Tij has suggested that stronger shell effects at N = 126 
and comparatively weaker shell effects at N — 82 are 
conducive to reproducing the r-process abundances in 
the second and the third peak, respectively, in the solar- 
system r-process abundance curve. Our results exhibit 
features which are consistent with this analysis and sup- 
port the conjecture that a different nature of the shell ef- 
fects may be at play in r-process nucleosynthesis of heavy 
nuclei. 

A scrutiny of the shell effects with various RMF forces 
and with various mass formulae has shown that a given 
nature of the shell effects in a known region may not 
be sufficient to calibrate the shell structure in unknown 
regions. Accordingly, an extrapolation in the unknown 
regions of the r-process path is frought with an uncer- 
tainty. 
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